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ABSTRACT 


A  completely  integrable  partial  differential  equation  is  one  which  has  a 
Lax  representation,  or,  more  precisely,  can  be  solved  via  a  linear  integral 
equation  of  Gel 1 f and-Levitan  type,  the  classic  example  being  the  Korteweg- 
de  Vries  equation.  An  ordinary  differential  equation  is  of  Painleve  type  if 
the  only  singularities  of  its  solutions  in  the  complex  plane  are  poles.  It  is 
shown  that,  under  certain  restrictions,  if  G  is  an  analytic,  regular 
symmetry  group  of  a  completely  integrable  partial  differential  equation,  then 
the  reduced  ordinary  differential  equation  for  the  G-invariant  solutions  is 
necessarily  of  Painleve  type.  This  gives  a  useful  necessary  condition  for 
complete  integrability ,  which  is  applied  to  investigate  the  integrability  of 
certain  generalizations  of  the  Korteweg-de  Vries  equation,  Klein-Gordon 
equations,  some  model  nonlinear  wave  equations  of  Whitham  and  Benjamin,  and 
the  BBM  equation. 
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SIGNIFICANCE  AND  EXPLANATION 


Several  very  important  nonlinear  partial  differential  equations  can  be 
solved  by  the  method  known  as  "inverse  scattering".  This  in  effect  reduces 
the  solution  of  the  nonlinear  equation  to  that  of  a  linear  integral  equation, 
and  it  is  obviously  desirable ,  given  any  nonlinear  partial  differential 
equation,  to  determine  whether  its  solution  is  amenable  to  this  technique. 
Hitherto  this  determination  has  been  largely  a  matter  of  chance,  but  the 
present  paper  gives  a  relatively  simple  systematic  test.  One  looks  at  the 
ordinary  differential  equations  satisfied  by  similarity  solutions  of  the 
nonlinear  equation.  If  these  ordinary  differential  equations  are  not  "of 
Painleve  type",  i.e.  if  they  possess  solutions  having  singularities  other 
than  poles,  then  the  non linear  equation  is  not  soluble  by  inverse  scattering. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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THE  CONNECTION  BETWEEN  PARTIAL  DIFFERENTIAL  EQUATIONS  SOLUBLE  BY  INVERSE 
SCATTERING  AND  ORDINARY  DIFFERENTIAL  EQUATIONS  OF  PAINLEVE  TYPE 
J.  B.  McLeod  and  P.  J.  Olver 


1 •  Introduction 

The  solution  of  certain  nonlinear  partial  differential  eouations  by  inverse  scattering 
techniques  has  been  the  fubject  of  considerable  interest  in  recent  years*  This  technique 
dates  back  to  a  fundamental  observation  of  P.  Lax,  [26] ,  that  some  puzzling  earlier  results 
of  Miura,  Gardner  and  Kruskal,  [33]  ,  on  the  relationship  between  the  Kortewecr-de  Vries 
equation  and  the  eigenvalue  problem  for  Hill’s  eouation,  could  be  placed  in  ar.  extremely 
lucid  and  general  theoretical  framework.  Lax’s  basic  idea  was  that  if  a  partial 
differential  equation  could  be  cast  into  the  form 

~~  =  [B,L]  *  BL  -  LB  ,  (1.1 ) 

(It 

where  L  and  R  are  linear  (differential)  operators  on  some  Hilbert  space  with  P  skew- 
adjoint,  then  the  eigenvalues  of  L  would  be  independent  of  time  t.  We  will  call  (1.1)  a 
Lax  representation  of  the  given  partial  differential  equation,  and  a  system  of  partial 
differential  equations  which  can  be  so  represented  is  often  called  completely  integrable, 
although,  as  will  be  seen  later  in  this  introduction,  we  shall  in  this  paper  find  it 
convenient  to  use  a  somewhat  different  definition  of  complete  integrability •  (The 
terminology  stems  from  the  interpretation  of  the  KdV  equation  as  a  completely  intearable 
Hamiltonian  system,  as  discovered  by  Gardner,  [16] ,  and  developed  in  great  detail  hv 
McKean,  van  Moerbeke  and  Trubowitz,  [30),  [31].) 

For  the  KdV  equation,  which  is 

u  +  6uu  +  u  =  0  ,  (1.2) 

t  X  XXX 


the  operators  appearing  in  the  Lax  representation  are 


w 


2 

L  ■  -D  -  u  , 

B  *  -{4D3  +  3 (Du  +  UD>}  , 

where  D  =  d/dx.  The  problem  of  finding  the  potential  u  from  the  spectral  data  of  L, 
known  as  the  "inverse  scattering  problem" ,  was  solved  by  Gel'fand  and  Levitan#  [19].  (See 
also  [14]#  [42]  for  comprehensive  introductions.)  Their  method  involves  the  solution  of  a 
certain  linear  integral  equation,  the  Gel 1 fand-Levitan  equation,  which  in  general  is  of  the 
following  form: 


K(x,y)  +  F(x,y)  +  /  K ( x , z ) H ( z ,y ) dz  *  0  .  (1.3) 

x 

(Technically,  (1.3)  is  the  Marchenko  form,  [29],  of  the  Gel 9 fand-Levitan  equation.)  Here 
F  and  H  are  constructed  from  the  relevant  spectral  data  of  L.  Once  K  has  been  found, 
the  potential  u  is  recovered  from  the  values  of  K  on  the  diagonal  x  *  y.  (In  the 
KdV  case,  u  =  2D (K(x,x) ] . )  Thus  the  basic  technique  for  solving  the  initial  value 
problem  for  a  completely  integrable  system  of  partial  differential  equations  consists  of 
the  following  steps: 

1)  Given  the  initial  data  u(x,0),  determine  the  appropriate  spectral  data  of  the 
operator  L  at  time  t  *  0. 

2)  Find  the  time  evolution  of  the  spectral  data,  and  hence  of  the  kernel  functions 
F  and  H  used  in  the  Gel 1 f and-Levitan  equation. 

3)  Solve  the  Gel ' fand-Levitan  equation,  regarding  t  as  a  parameter,  and  thus 
recover  the  solution  u  of  the  oricfinal  system. 

The  inverse  ncatterina  technique  outlined  above  has  been  applied  to  a  number  of 
physically  relevant  partial  differential  equations.  Of  particular  interest  is  the  work  of 
Sakharov,  ManaKov  and  Shabat,  [48] -[53],  on  the  nonlinear  Schrodinger  and  other  physically 
interesting  equations,  and  of  Ablowitz,  Kaun,  Newell  and  Segur  on  a  general  class  of 
2x2  matrix  systems,  [2],  and  the  three-wave  interaction  equations,  [24].  In  all  cases 
the  appropriate  Gel ' fand-Levitan  equation  assumes  the  form  (1.3),  although  F,  H  and  K 
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may  be  matrix- valued  functions.  The  formula  for  recovering  the  solution  u  from  the 


kernel  K  varies  from  case  to  case. 

The  one  notable  drawback  in  this  elaborate  theory  is  that  for  a  given  system  of 
partial  differential  equations  there  has  not  to  date  been  any  systematic  method  of 
determining  whether  or  not  it  has  a  Lax  representation  and,  if  so,  how  to  find  the 
operators  B  and  L.  Previous  work  has  relied  either  on  inspired  guesswork,  or  else  on 
fixing  the  form  of  the  operators  B  and  L  and  seeing  what  systems  of  partial 
differential  equations  result.  In  particular,  recent  work  of  Gel'fand  and  Dikii,  [17]  # 
[18],  which  lists  several  other  references,  has  used  abstract  differential-algebraic 
methods  to  give  a  general  classification  of  pairs  of  differential  operators  B  and  L 
such  that  (1.1)  corresponds  to  a  bona  fide  system  of  partial  differential  equations.  This 
approach,  however,  while  providing  a  large  number  of  completely  integrable  systems,  is  not 
yet  able  to  answer  the  above  question  of  whether  a  given  system  is  completely  integrable. 

The  basis  of  the  present  paper  is  an  observation  of  Ablowitz  and  Segur,  [5] ,  that  the 
equations  for  the  group- invariant  ( self- similar )  solutions  of  known  examples  of  completely 
integrable  equations  turn  out  to  be  ordinary  differential  equations  studied  extensively  by 
Painleve  and  his  students  around  the  turn  of  the  century,  [10],  [38].  (See  also  [22], 

[23]  for  general  accounts  of  the  subject  and  further  references.)  These  equations  are 
characterized  by  the  property  that  all  their  solutions  in  the  finite  complex  plane  possess 
only  poles  as  singularities,  and  hereafter  we  will  refer  to  an  ordinary  differential 
equation  with  this  property  as  an  equation  of  Painleve  type.  (Painleve  allowed  also  fixed 
singularities  of  an  arbitrary  type,  but  we  will  not.)  Hastings  and  McLeod,  [21],  exploited 
this  relationship  to  solve  a  nonlinear  connection  problem  for  the  second  Painleve 
transcendent,  and  conjectured  that  the  above  relationship  was  not  fortuitous: 

Conjecture .  If  a  system  of  partial  differential  equations  is  completely  integrable, 
and  G  is  a  symmetry  group  of  this  system,  then  the  reduced  system  of  ordinary 
differential  equations  for  the  G-invariant  solutions  is  of  Painleve  type. 

This  conjecture,  if  true,  would  provide  a  powerful  necessary  condition  to  test  for 
complete  integrability .  Here  we  will  prove  a  somewhat  weakened  version  of  the  conjecture. 


which  nevertheless  proves  useful  in  several  applications.  There  are  two  restrictions. 

First,  if,  in  the  Lax  operator  L,  some  combination  of  the  solution  u  and  its  spatial 

derivatives  occurs,  say  Q(u),  then  it  is  this  combination  (or  combinations)  that  must 

2 

have  only  poles  as  singularities.  For  instance,  if  L  =  D  +  ux »  then  only  u  is 
required  to  have  poles,  and  thus  we  may  allow  logarithmic  branch  points  as  singularities  of 
the  solutions  of  the  reduced  ordinary  differential  equations.  Usually  we  will  assume 
that  Q  is  a  linear  combination  of  u  and  its  spatial  derivatives,  callinq  this  case 
linearly  completely  infegrable.  Secondly,  the  same  combination  Q  must  satisfy  certain 
preconditions  for  the  inverse  scattering  formalism  to  go  through?  this  means  that,  when 
restricted  to  the  real  axis,  Q  either  is  periodic  or  satisfies  decay  conditions  at 
x  -  +  09 ,  which  implies  corresponding  restrictions  on  the  solutions  u  that  can  be 
considered.  It  is  only  for  such  solutions  that  Q(u)  must  be  meromorphic.  If  a  system  of 
ordinary  differential  equations  has  the  property  that,  for  such  solutions  u,  the 
combination  Q(u)  is  meromorphic,  we  say  that  the  system  is  of  restricted  Painleve  type 
relative  to  Q.  Our  basic  result,  in  rough  form,  replaces  "Painleve  type"  by 
"restricted  Painleve  type"  in  the  above  conjecture. 

The  first  requirement  for  stating  and  proving  a  precise  form  of  the  conjecture  is  to 
define  what  is  meant  by  a  partial  differential  equation  beina  completely  integrahle. 

Father  than  take  the  Lax  representation  as  our  starting  point,  we  shall  adopt  the  more 
practical  view  of  the  Gel  * fand-Levitan  equation  being  of  primary  importance.  Thus  a 
completely  integrable  system  of  partial  differential  equations  is  defined  as  one  w'-iose 
solutions  are  found  by  solvinq  a  linear  integral  equation  of  special  type,  cp.  Definition 
2.1  below.  This  viewpoint  is  necessitated  by  the  fact  that  the  formal  Lax  representation 
theory  of  Gel’fand  and  Dikii  has  not  reached  the  point  of  statinq  an  analoque  of  the  result 
of  Gel'fand  ami  Levitan  for  an  arbitrary  scatterinq  operator  L,  although  some  formal 
progress  in  this  direction  has  been  made  in  [28]  ,  [53].  Since  in  all  examples  known  to  th*-* 
authors  the  Gel ' fand-Levitan  equation  is  always  of  the  form  (T.3),  it  might  seem  reasonable 
that  this  should  be  true  in  general,  but  Application  V  in  Section  4  below  indicates  that  it 
may  be  necessary  to  restrict  L  at  least  to  being  of  prime  order. 


ts  tt0S? Sf9 


The  main  tool  in  our  proof  is  a  theorem  of  Steinberg,  [43],  which  states  that  if 
T(z)  is  an  analytic  family  of  compact  operators  in  a  Banach  space,  then  (I  -  T(z))  1 , 
provided  this  inverse  exists  for  at  least  one  value  of  z,  is  a  meromorphic  family  of 
operators.  Under  appropriate  assumptions  on  the  initial  data  of  our  completely  integrable 
system  to  ensure  that  the  functions  F  and  H  in  the  Gel ’ f and-Levitan  equation  satisfy 
certain  analyticity  criteria,  we  can  conclude  from  Steinberg's  result  that  Q  must  be  a 
meromorphic  function  of  (x,t).  Now  suppose  that  G  is  a  one-parameter,  analytic,  regular 
local  group  of  transformations  acting  on  the  space  of  independent  and  dependent  variables 
which  leaves  the  set  of  solutions  of  the  system  of  partial  differential  equations 
invariant.  Then  the  G-invariant  (self-similar)  solutions  can  all  be  found  by  integrating  a 
system  of  ordinary  differential  equations  on  the  quotient  manifold  whose  points  correspond 
to  the  orbits  of  G.  (See  Section  3  for  details.)  The  analyticity  of  G  implies  that  for 
any  G-invariant  solution  whose  initial  data  satisfies  the  inverse  scattering  assumptions, 
the  function  Q  on  the  quotient  manifold  can  have  only  poles  for  singularities.  In  other 
words,  the  reduced  system  of  ordinary  differential  equations  must  be  of 

restricted  Painleve  type  relative  to  p.  This  completes  the  outline  of  the  proof  of  our 
main  theorem;  precise  statements  and  proofs  will  appear  in  Sections  2  and  3. 

In  Section  4  we  discuss  some  applications  of  this  result.  First  we  show  that  the 
generalized  KdV  equation 

u  +  uPu  +  u  =0 

t  X  XXX 

can  be  linearly  completely  inteqrable  only  if  p  =  0,  1,  or  2.  These  exceptional  cases 
correspond  to  the  Airy  equation  in  moving  coordinates,  the  KdV  and  the  modified  KdV 
equations,  which  are  well  known  to  be  completely  integrable.  Secondly  we  consider  a 
nonlinear  Klein-Gordon  equation  in  characteristic  coordinates 

u  =  f 1 ( u )  .  (1-4? 

xt 

It  is  shown  that  if  f(u)  is  a  rational  function,  real  for  real  u  and  with  two 
consecutive  zeros,  simple  or  double,  on  the  real  axis,  and  if  (1.4)  is  linearly  completely 
integrable,  then  f  is  a  polynomial  of  degree  at  most  4.  Further,  if  f(u)  is  a  linear 
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a  .u 

combination  of  exponentials  e  3  with  the  all  rational  multiples  of  some  complex 

number  a,  again  real  for  real  u  and  with  two  consecutive  simple  or  double  zeros,  and  if 

(1.4)  is  linearly  completely  integrable,  then 

r/  .  2  0u  0u  .  -Bu  ^  -20u 

f(u)  =  c2e  +  c^e  +  cQ  +  c^e  +  c_2e  * 

for  some  number  0.  The  next  application  shows  that  certain  nonlinear  model  wave  equations 

considered  by  Benjamin,  Bona  and  Mahony,  [7],  and  Whitham,  [46],  cannot  be  linearly 

completely  integrable.  The  last  example  deals  with  the  BBM  equation,  [7] , 

u.  +  uu  -  u  .  *  0  .  (1.5) 

t  x  xxt 

Although  this  cannot  be  treated  rigorously  by  the  methods  of  the  present  paper,  we  show 
that  if  the  full  conjecture  were  true,  then  (1.5)  could  not  be  linearly  completely 
integrable.  From  these  results,  it  can  be  seen  that  our  criterion  for  complete 
integrability  is  a  powerful  preliminary  test  when  considering  whether  or  not  a  system  of 
partial  differential  equations  can  be  integrated  by  inverse  scattering  techniques. 

A  recent  preprint  of  Jblowitz,  Ramani  and  Segur,  [3],  also  considers  the  above 
conjecture.  They  prove  a  similar  result,  although  they  place  a  rather  severe  restriction 
on  the  form  of  the  functions  occurring  in  the  Gel' fand-Levitan  equation  for  the  group- 
invariant  solutions  which  is  not  sufficiently  justified.  Also,  the  only  groups  they 
consider  are  groups  of  scale  transformations;  the  more  general  groups  that  we  consider 


allow  a  much  wider  range  of  applications 


2. 
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Analyticity  Properties  of  Completely  Integrable  Differential  ^uations 
Consider  a  system  of  partial  differential  eauations 

A( t /X / u)  =  n  ,  (2.1; 

where  x,t  £  R  and  u  -  (u\...,um)  £  Rm  is  a  vector-valued  function*  We  assume  that  the 
initial  value  problem  of  (2.1)  with 

u(x,0)  =  f ( x )  (2.2) 

is  well  posed  for  f  in  some  Banach  space  8  of  functions,  so  that  for  t  sufficiently 
small,  there  is  a  unique  solution  u(x,t)  of  (2.1-2).  In  practice  8  is  either  a  space 
of  functions  decreasing  sufficiently  rapidly  at  +  «  or  a  space  of  periodic  functions. 
Usually  the  presence  of  appropriate  conservation  laws  will  ensure  that  the  solutions  are 
actually  global  in  t,  but  this  will  not  be  assumed  a  priori*  The  first  task  is  to  make 
precise  what  is  meant  by  (2.1)  being  completely  integrable.  Rather  than  use  the  usual  Lax 
representa tion  of  the  ecmation,  we  will  assume  a  more  practical  outlook  and  take  the 
Gel' fand-Levitan  integral  equation  as  our  starting  point.  As  the  examples  demonstrate,  in 
all  known  cases  there  is  such  an  integral  ecruation  for  solution  of  the  inverse  scattering 
problem  for  the  operator  L  in  the  Lax  representation  (1.1). 

Definition  2.1.  A  system  of  partial  differential  eauations  is  completely  integrable 
relative  to  £(u)  in  the  Panach  space  8  if  there  is  a  linear  matrix  integral  equation  of 
the  form 


00 

K( x ,y ;t )  +  F( x , v ; t )  +  /  K ( x , z ; t ) H ( z , y ;t ) dz  =  0  ,  (2.1) 

x 


called  the  Gel* fand-Levitan  equation,  satisfying  the  following  properties: 

i)  F,  H,  V  are  fl  x  tj  matrices  of  functions; 

ii)  r  and  H  are  uniquely  determined  by  the  initial  data  (2.2); 

iii)  for  initial  data  in  8,  and  for  all  real  x,  y,  all  complex  c,  and  t  in 
some  domain  7i  in  C,  the  functions  F(  x  -  et,v  -  et;t)  and  H(x  -  et,y  -  ct;t)  are 

analytic  in  £,  t,  and  there  is  a  Banach  space  8  (not  necessarily  the  same  as  F)  for 

* 

which  Ff  x  -  ct  ,v  -  ft  j  t )  £  8  as  a  ^unction  of  y  and  the  operator 
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et,y  -  et;t)dz 


T(x,t)f(y)  =  /  f ( z)H( z  - 
x 

_  * 

is  a  compact  operator  in  o  ; 

iv)  the  Gel ' fand-Levitan  equation  has  a  unique  solution  (in  8*)  for  all  x  and  at 
least  one  t  in  ft; 

v)  the  solution  u  of  the  system  (2.1-2)  can  he  recovered  from  the  solution  K  of 
the  Gel' fand-Levitan  equation  via  a  relation  of  the  form 

Q[u(x,t)]  =  P[K(x,x,t)] ,  (2.4) 

where  Q  is  some  function  of  u  and  its  spatial  derivatives  and  P  is  a  polynomial  in 
y  and  its  spatial  derivatives. 

Thus  to  recover  the  solution  u  of  a  completely  integrable  system  of  partial 
differential  equations,  we  must  solve  the  Gel' fand-Levitan  equation  for  K,  and  then  solve 
the  differential  eouation  (2.4)  for  u.  In  practical  examples,  Q  is  a  linear  combination 
of  the  spatial  derivatives  of  u,  and  in  this  case  the  system  will  he  called  linearly 
completely  integrable.  It  should  also  be  remarked  that  the  requirement  that  iii)  hold  for 
all  complex  e  can  certainly  be  relaxed,  although  there  seems  little  practical  point  in 
doing  90,  and  that  the  domain  ft  will  customarily  include  the  origin  or  at  least  have  the 
origin  on  its  boundary  (it  might,  as  in  the  example  of  the  Kdv  eauation  below,  be  a  sector 
of  a  circle  centre  the  origin) • 

We  now  illustrate  the  definition  with  two  well  known  examples. 

Example  2.2:  The  Korteweg-de  Vries  Equation 

This  is  the  original  example  of  the  use  of  inverse  scattering  techniques,  [2f ) , 

[32].  The  equation  is 

u  +  f»uu  +  u  *  0  ,  ( 2  •  5 ) 

t  X  XXX 

and  has  a  Lax  representation  with  operators 

L  -  -D2  -  u  ,  (2.*) 

B  »  -{4D3  +  3 (Du  +  UP) }  , 
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where  D  =  d/dx.  In  the  case  that  these  operators  act  on  a  dense  subspace  of  l/(R)  the 

appropriate  spectral  data  of  L  consists  of 

2  2 

i)  the  eiaenvalues  A.=-k A  =  -k  , 
li  n  n 

ii)  the  associated  norming  constants,  c-,...,^,  defined  so  that  if  is  the 

-V  .  x  J 

eigenfunction  associated  with  A.  satisfying  «£.(x)  ~  e  **  ,  x  +®,  then 
00  D  1 

-1  r  .2, 

c  .  -  J  v?  .  dx  , 

3  :»  ^ 

2 

iii)  the  reflection  coefficient  R(k)  for  k  real  and  positive,  k  *  A,  defined 
so  that  there  exists  a  solution  y  of  Ls?  =  Ay  which  satisfies  the  boundary  conditions 


-*ikx  ,  ikx 

e  +  *>{k)e  , 


x  +  +«  , 


-  ikx 


'  ( x)  - 

T(k)  e 

The  time  evolution  of  the  spectral  data  is  determined  by 

A .  ( t)  =  X.(0)  , 

2  2 

c.(t)  =  c.(0)  expf8k^t]  , 

2  2  3 

R(k,t)  =  R(k,0)  exp[Rik^t]  * 


f  2.7) 


The  Gel ' f and-Levi tan  equation  takes  the  form 


k(x,y;t)  +  F( x  +  y;t)  +  /  X(x,z;t)F(z  +  y,t)dz  =  0  , 


(?.*) 


and  the  solution  of  the  KdV  equation  is  recovered  via  the  simple  formula 

u(x,t)  -  2  ~  K(x,x;t)  • 
dx 

The  kernel  F  of  the  Gel'  fand-Levitan  equation  is  criven  by 


(  ?  .  «  ) 


F(  x  ,  t )  =  F  c  .  ( t)  exp  [-k  .  x]  +  /  R(k,t)e^Xdk  ,  <  2 .  1  n' ' 

1  2  2ti 


j  =  1 


and  this  solution  is  val  id  provided  the  initial  data  u(x,0)  =  f(  x)  satisfies 


[  (1  +  |x|)|f(x)|dx  <  00  . 


(2.11) 


_o  _ 


The  uniqueness  of  the  solution  of  (2.8)  in  the  KdV  case  is  a  standard  result ,  and  the 
only  item  remaining  to  be  checked  is  Definition  2.1  in  condition  iii).  So  far  as 
analyticity  is  concerned,  the  only  part  of  F  that  could  fail  to  be  analytic  is  that 
corresponding  to  the  continuous  spectrum  of  L: 


F  (x,t) 
c 


/  R(k,t)  eikXdk 
-00 

00 

/  R(k,0 )exp[8ik3t  +  ikx] dk  . 


If  we  take  any  reasonable  space  of  initial  data  for  8,  for  example  that  given  by  (2.11), 

then  R(k,0)  can  be  extended  analytically  into  the  upper  half  of  the  k-plane,  and 
2 

(R(k,0)/k  !  is  bounded  as  |k!  ♦  °®.  (The  function  R  is  closely  related  to  the  spectral 
density  function  m  of  Titchmarsh,  and  the  analyticity  and  estimates  can  be  obtained  by 
suitably  translating  the  results  in  [44,  Chapter  V).)  If  therefore  we  write 

0  * 

( x , t )  *  —  {/  +  /  }p(k, 0 }exp[8ik  t  +  ikx] dk  *  F^  +  F^ ,  say  , 

-00  0  ^ 

and  consider  sufficiently  F^ ,  then,  if  t  is  real  and  positive,  we  can  deform,  the 
integral  from  (O,00)  to  (0,a°°),  for  any  a  with  0  <  arg  a  <  —  tt.  We  can  now 
increase  arg  t,  but  the  range  for  a  becomes  0  <  arg  a  <  —  (tt  -  arg  t).  Nonetheless 
this  does  allow  us  to  define  Fj(x,t)  as  an  analytic  function  of  t  for  0  <  arg  t  <  tt. 
(It  is  also  an  analytic  function  of  x  since  for  large  k  the  term  k^t  dominates 
kx.)  If  we  decrease  arg  t,  the  range  for  a  becomes  -  —  arg  t  <  arg  a  <  -j  tt,  which 
allows  us  to  define  F,|(x,t)  as  an  analytic  function  of  t  for  -tt  <  arg  t  <  0,  and  so 
in  fact  in  the  whole  complex  plane  cut  along  the  negative  axis.  Similar  remarks  apply  to 

v 

Further,  by  using  the  deformed  contours  and  integrating  by  parts  (integrating  e 
and  differentiating  the  remainder),  we  see  that  F  e  8,  the  Banach  space  defined  by 
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(2.11),  and  that  the  operator  T  is  compact  in  R,  although  8  is  certainly  not  the  only 
possible  choice  for  8*. 


Example  2.3:  The  AKNS  Systems 

A  generalized  eigenvalue  problem  considered  by  Ablowitz,  Kaup,  Newell  and  Segur ,  [1]# 
[21  ,  continuing  work  of  Zakharov  and  Shabat,  [51] ,  assumes  the  following  form: 


v^  *  -i£v  ♦  q(x)w  , 

w  *  i£w  +  r(x)v  , 
x 


(2.12) 


where  £  corresponds  to  the  eicrenvalue  and  q{ x) , r ( x)  +  0  as  |x|  ♦  00 .  If  we  include  a 
"time"  dependence  and  suppose  that  the  potentials  g(x,t),  r(x,t)  evolve  according  to  some 
evolution  equation  in  such  a  way  that  £  is  independent  of  t  and,  say, 


v  »  a(x,t,£)v  +  B(  x, t ,  £)w  , 
w  *  Y (x , t , £ ) v  -  a(x,t,£)w  , 


(2.13) 


then  cross-differentiation  between  (2.12)  and  (2.13)  leads  to 

a  *  -rB  +  ay  , 
x 

3x  +  2i £ 8  =  q  -  2qa  ,  (2.14) 

-  2i£y  *  r  +  2ra  . 

The  equations  (2.14)  can  also  be  reached  via  a  Lax  representation  (1.1),  if  we  take 


r  iD 

-iq 

,  *  * 

'a  6  ] 

i 

ir 

-iD  j 

I 

Y  -a  J 

The  insistence  that  (1.1)  should  be  valid  when  operating  on  solutions  of  (2.12)  leads 

again  to  (2.14),  and  conversely. 

Ry  judicious  choice  of  the  functions  a,  8,  Y,  the  AKNS  system  (2.14)  gives  rise  tc  a 
larae  variety  of  physically  interesting  equations,  including  the  KdV,  modified  KdV,  sine- 
Gordon,  nonlinear  Schrodinger,  and  others.  For  instance, 


q  «  U,  r  *  u  , 

a  -  -2U2  -  i  |u |  2  , 
8  -  2Cu  +  iu^  , 

r  =  2c u  -  iux  , 


leads  to  the  nonlinear  Schrodinger  equation 

2 

iu  =  -u  +  2u|u|  , 

t  XX 

whereas 


q 

a 

8 

leads  to  the  sine-Gordon  equation 


2  Ux' 


<c  —  -  ' 

Y  “  4C  8in 


u  “  sin  u  ♦ 
xt 

In  fact,  a  system  of  AKNS  type  with  any  given  dispersion  relation  can  be  found! 

The  spectral  data  for  (2.12)  is  constructed  as  follows.  For  C  real,  consider  the 
solutions  satisfying  the  following  asymptotic  boundary  conditions: 


* 

* 


-iCx 

** 

r  °i  U* 

v?  - 

Ux 

■  t 

x  ♦  -®  ; 

x  ♦  . 


Then 

*  -  a(Z,t)ty  +  b(c,t>^  , 

*  «  -a( C ,t ) ^  +  b(£,t)^ 


for  some  a,b,a,b.  If  (1  +  |x|)q,(1  +  Jx|)r  e  l\r',  then  a  ( respectively  a)  can,  as 
functions  of  5,  be  analytically  extended  to  the  upper  (lower)  half  plane,  and  their  zeros 
correspond  to  the  eigenvalues  of  (2.12).  Let  F  =  Fc  ♦  F^,  where 


F  ( x ,  t ) 
c 


i_  r 

2it  1  a(C,t) 


eliXdC 
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which  corresponds  to  the  continuous  spectrum,  and 


N  i;  x 

F  (x,t)  *  -i  y  c  e  n 
d  n 

n*1 


which  corresponds  to  the  discrete  spectrum  consisting  of  the  zeros. 


of  a(C)  in  the 


upper  half  plane  (for  convenience  assumed  simple)  and  appropriate  norming  constants  c  ^ . 

Similarly  construct  F  =  F  +  F,  from  a  and  b.  Then  if 

c  d 


-F 

0 


the  Gel' f and-Levi tan  equation  takes  the  form  (with  t-dependence  suppressed) 


K(x,y)  +  F(x  +  y )  +  /  K(x,z)F(z  +  y)dz  =  0  , 
x 


where 


The  potentials  q  and  r  are  recovered  via 

q(x)  =  -2K1(x,x),  r(x)  =  -2K2(x,x)  , 

(Note  that  the  sine-Gordon  case  gives  an  example  in  which  the  operator  Q  in  (2*4)  is 

nontrivial.)  Finally,  it  remains  to  discuss  the  time  evolution  of  the  spectral  data. 

Under  the  assumption  that 

a(x,t,C)  ♦  a  (C),  6(x,t,C)  ♦  0,  *Y(x,t,C)  ♦  0 

o 

as  |  x  |  +  *»,  it  can  be  shown  that 


aU,t)  *  aU,0)  , 

b(  C / 1 )  =  b(C,0)exp[-2ao(c)tl  , 

C^t)  -  C.(0)  , 

cJ(t)  *  c,(0)exp(-2a  (C.)t)  • 
1  j  o  j 


ny  . . 


If,  in  addition,  I  Ot^  (  ^ )  /  C I  +  00  as  I  C I  *  00  /  then  arguments  similar  to  those  used  in  the 
previous  example  show  that,  for  suitable  initial  data,  f(x,t)  is  analvtic  in  t  for 
all  x,  and  also  analytic  in  x.  a  larae  class  of  AKNS  svstems  satisfies  these  additional 
assumptions . 

We  now  investigate  the  properties  of  the  solutions  of  a  general  intearal  equation  of 
Gel ' f and-Levitan  type.  Our  main  tool  is  the  following  theorem  of  Steinberg,  [43], 
general izincr  a  theorem  of  Dolph,  McLeod  and  Thoe,  [13],  for  the  case  of  Hi lber t- Schmidt 
operators . 

Theorem  2.4.  Let  8  be  a  Banach  space,  and  let  T(z)  be  an  analytic  family  of 
compact  operators  defined  for  z  e  £]  C  C.  Then  either  I  -  T(z)  is  nowhere  invertible  for 

z  £  ij  or  (I  -  T(z))”1  is  meromorphic  for  z  e  ft. 

Let  us  write  the  Gel ' f and-Levitan  enuation  (2.3)  in  the  symbolic  form 

(I  +  T ( x ? t ) )K(x,y;t)  +  F(x,y;t)  =  0  ,  (2.15) 

where  T(x,t)  denotes  the  family  of  intearal  operators 

00 

T( x ,t ) f { y )  =  /  f( z)H{z,v;t)dz  .  (2.16) 

x 

It  will  always  be  assumed  that  T(x,t)  is  a  compact  operator  for  each  fixed  (x,t).  For 
instance,  this  is  guaranteed  if 


00  oo 

/  /  lH(v,z?t)  |  2dyd2  <  «>  , 

X  X 

indeed,  in  this  case  T  is  Hilbert-Schmidt . 

To  apoly  Steinberg's  theorem,  we  treat  the  time  t  as  the  complex  parameter.  (Note 
that  it  would  not  do  any  aood  to  look  at  x  as  this  parameter  since  the  domain  of 
integration  for  T(x,t)  depends  on  x,  and  so  the  operators  could  not  possibly  be 
analytic  for  a  large  enough  class  of  functions.)  Now,  for  all  x,  y,  if  the  kernel 
H(x,y?t)  depends  analytically  on  t  for  t  e  then  the  operators  T(x,t)  depend 

analytically  on  t.  Tf  furthermore  F(x,y;t)  is  analytic  in  t,  then  Steinberg's  theorem 
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implies  that 


i 


K(x,y?t)  =  -{I  +  T(x,t))”1F(x,y;t) 

is#  for  each  fixed  (x,y),  a  meromorphic  function  of  t.  (It  is  one  of  the  assumptions  of 
complete  integrabil ity  that  the  inverse  exists  for  at  least  one  t.)  Therefore 

Q [u{ x , t )  ]  =  P[K(x,x;t)J 

is  also  a  meromorphic  function  of  t  for  each  fixed  x. 

Theorem  2*5.  If  a  system  of  partial  differential  equations  is  Q- completely  integrable 
in  the  Banach  space  8,  and  if  the  initial  data  u(x,0)  e  8#  then  t-he  function 
Q [ u( x , t ) ]  is  meromorphic  in  t  for  t  e  ft  and  each  fixed  x. 

A  slight  generalization  of  this  theorem  will  prove  to  be  of  use  in  the  sequel* 

Suppose  that  the  time  axis  is  "skewed",  by  making  the  change  of  variables 

(x#t)  ■  (x  +  et,t) 

for  some  real  e.  If  u  =  f(x,t)  is  the  solution  to  the  "unskewed"  equation#  then 
u  *  f(x,t)  =  f(x  -  £ t# t )  is  the  solution  in  terms  of  the  new  coordinates.  If  we  let 

K(x#y;t)  =  K(x  -  et,y  -  £t;t)  # 

then  K  is  a  solution  of  a  Gel' fand-Levitan  equation  of  the  form 


K(x,y?t)  +  F(x  -  Et#y  -  et;t)  +  /  K(x,z;t)H(z 

x 


et#y  -  et; t ) dz  -  0 


Therefore  the  "skewed"  equation  is  also  completely  integrable#  which  gives  the  following 
theorem. 

Theorem  2.6.  If  a  system  of  partial  differential  equations  is  Q-completely  integrable 
in  the  Banach  space  8,  and  if  the  initial  data  are  in  8#  then  the  function  Q[u(x,t)] 
is  meromorphic  in  (x#t)  for  x  r.  C,  t  e  fi. 
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3 .  Symmetry  Groups  and  Group-Invariant  Solutions 

Given  a  system  of  partial  differential  equations 

A  { x  #  u )  =  0,  x  e  Rm,  u  e  Rn  #  (3.1) 

a  symmetry  group  will  be  a  local  Lie  group  of  transformations  G  acting  on  the  space 

Rm  x  Rn  of  independent  and  dependent  variables  which  preserves  the  set  of  solutions  of 

(3.1).  The  group  acts  on  the  solutions  by  transforming  their  graphs.  For  simplicity#  we 
will  restrict  our  attention  to  projectable  groups#  meaning  those  in  which  the 
transformations  are  all  of  the  form  (x#u)  =  (a(x) #$(x#u) ) .  If  u  *  f(x)  is  a  solution, 
the  transformed  solution  will  be  given  by  the  formula 

u  *  f(x)  =  B(a  ^  ( x )  #  f  ( a  Vx)))  , 

provided  a  is  invertible.  The  details  of  this  theory  can  be  found  in  [8] ,  [34]  and  [37]  . 

If  G  is  a  local#  projectable#  one-parameter  group#  its  infinitesimal  generator  is  a 
vector  field  of  the  form 


hi  ,  #v  n 

ri.  .  3  .  t  1# 


v  *  l  r<x>  +  l  ^<x#u)  — 

i«=l  3x  j=1  3uJ 


Denote  the  partial  derivatives  of  the  uJ  by 


.  ,\y\  j 

3  « _ u 


“k  =  ,  k  k  ”  (ki  "••'V  ' 


ax" 


and  let 


3  r  j  3 

Di  *  T7  +  /  Uk#i  ' 

3x  j  #  k  3u£ 


where  k,i  -  (k1#...#k^  ♦  1,...,km),  denote  the  total  derivative  with  respect  to  x1. 
Define  the  prolongation  of  v  to  be  the  vector  field 


pr 


-*  “h  r  j  3 

V  -  V  4  l  *  l  —  # 


jr*  k  3UV 


(3.2) 


where 


2  m 

(Here  D  «  D  D  ...  D  ,)  Then  the  following  theorem  provides  the  infinitesimal 

i  2  m 

criterion  for  G  to  be  a  symmetry  group  of  a  given  partial  differential  equation. 

Theorem  3.1.  Let  G  be  a  connected  local  Lie  group.  Then  in  general  G  is  a 
symmetry  group  of  the  system  of  partial  differential  equations  (3.1)  if  and  only  if 

pr  v(  A)  -  0  (3.3) 

whenever  A  *  0  for  all  infinitesimal  generators  v  o£  G. 

A  precise  statement  and  proof  of  this  may  be  found  in  [34] ,  or,  in  simplified  form,  in 


[35]  . 


Example  3.2.  Consider  the  generalized  KdV  equation 

A(x,u)  Hu+u^u  +  u  *0, 

t  X  XXX 

where  p  is  a  positive  integer*  Let  G  denote  the  one-parameter  group  of  scale 
transformations 


(x,t,u)  ♦  (e  P^x,e  3p\,e2*u)  ,  X  £  R  . 


The  infinitesimal  generator  of  G  is 


v  =  -px  - - 3pt  —  +  2u  r-  . 

dx  dt  du 

From  the  prolongation  formula  (3.2),  we  find  that 

— ^  — *■  3  3  3  3 

pr  v  *  v  +  (p+2)u  t —  +  (3p+2)u  -r —  +  (2p+2)u  r -  +  ( 3p+2 )  t - 

*  r  x  3u  r  t  3u  xx  du  du 

X  t  XX  XXX 


(3.4) 


Therefore 

pr  v ( A]  *  (3p+2)ut  +  (3p+2)uPux  +  (3p+2)uxxx  *  (3p+2)A  , 
so  that  the  inf initesimal  criterion  (3.3)  is  verified,  we  conclude  that  if  u  =  f(x,t)  is 
any  solution  of  (3.4),  so  also  is  u  =  e^f  ( eP*x,e3p^t )  .  This  may  be  checked  directly. 

In  general,  to  find  the  symmetry  group  of  a  given  system  of  partial  differential 
equations,  one  looks  for  all  vector  fields  v  satisfying  (3.3).  This  necessitates  the 
solution  of  a  large  number  of  elementary  partial  differential  equations  for  the  coefficient 
functions  of  v,  which  is  easily  done  in  practice,  as  in  the  above  mentioned 

references . 
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Now,  given  a  symmetry  group  G,  a  G-invariant  (or  self-similar)  solution  of  (3.1)  is 
a  solution  which  is  unchanged  by  the  transformations  in  G.  The  fundamental  property  of  G- 


invariant  solutions  is  that,  roughly  speaking,  they  may  all  be  found  via  the  integration  of 
a  system  of  partial  differential  equations  in  fewer  independent  variables.  To  make  this 
precise,  we  must  assume  that  G  acts  "regularly”  in  the  sense  of  Palais,  [39].  This 
requires 

i)  that  all  the  orbits  of  G  have  the  same  dimension,  I, 
ii)  that,  for  any  point  (x,u),  there  exist  arbitrarily  small  neighborhoods  N  such 
that  the  intersection  of  any  orbit  0  of  G  with  N  is  a  pathwise  connected  subset  of  0. 

(The  prototypical  group  actions  excluded  by  the  second  requirement  are  the  irrational 
flows  on  the  torus.) 

Under  these  two  assumptions,  it  is  well  known  that  the  quotient  space  M  *  Rm  x  Rn/G, 
whose  points  correspond  to  the  orbits  of  G,  can  be  naturally  endowed  with  the  structure 
of  a  smooth  (although  not  always  Hausdorff)  manifold.  Moreover  the  G-invariant  solutions 
of  (3.1)  are  all  obtained  by  integrating  a  reduced  system  A/G  *  0  of  partial  differential 
equations  on  M,  which  necessarily  has  \  fewer  independent  variables.  Precise 
statements  and  proofs  of  these  results  may  be  found  in  [34] • 

For  our  purposes,  the  construction  of  the  reduced  system  for  the  G-invariant  solutions 
proceeds  as  follows:  Local  coordinate  systems  on  the  quotient  manifold  M  are  provided  by 
a  "complete  set  of  functionally  independent  invariants  of  G",  cf.  [37].  if  G  is 
projectable,  these  are  functions  of  the  form 

s  ( x) , . . • , £  (x),  w  (x,u),...,w  ( X , u )  , 

which  are  unchanged  under  the  action  of  G.  The  functional  independence  means  that  the 


Jacobian  matrix 


is  everywhere  nonsingular.  The  reduced  system  A/G  =  0  will  then  be  found  in  terms  of  the 


new  independent  variables  F,1  and  the  new  dependent  variables  w-*. 

Example  3.3.  Return  to  the  generalized  KdV  equation  and  the  group  of  scale 
transformations  G  of  the  previous  example.  Note  that  G  acts  regularly  on 

-  R3  -  {0},  the  origin  being  a  singular  point.  The  quotient  manifold  M  *  ZQ/G  is  a 
complicated,  non-Hausdorf f  manifold.  However,  if  we  restrict  our  attention  to  the  subset 
defined  by  t  >  0,  then  a  complete  set  of  functionally  independent  invariants  of  G  is 
provided  by 


Viewing  w  as  a  function  of 


5  - 

Z ,  we  compute 


2/3p 
t  u 


(3.5) 


u 


u 


t 


u 

X 


t-2/3”w 


-1 


u  =  t 


XXX 

the  primes  denoting  derivative^  with  respect  to  £.  Substituting  these  expressions  into 

(3.4),  and  factoring  out  t  3p,  we  obtain  the  reduced  ordinary  differential  equation 

-  w  -  ~  £w*  +  w^w'  +  w*  ' '  =0  .  (3.6) 

3p  3 

Every  solution  of  (3.6)  gives  rise  to  a  scale-invariant  solution  of  the  generalized  KdV 

equation,  via  the  transformation  (3.5).  In  the  special  case  p  =  2  (the  modified  KdV 

equation),  (3.6)  can  be  integrated  once,  so  that 

.  .  1  3  1  r 

w”  =  -  j  w  +  —  £w  +  k 

for  some  constant  k.  This  may  be  recognized  as  the  second  Painleve  transcendent,  [23]. 

Now  we  restrict  our  attention  to  a  Q-completely  integrable  system,  A  =  0,  of  partial 
differential  equations  in  two  independent  variables,  (x,t).  Let  G  be  a  one-parameter 
local  projectable  symmetry  group  of  the  given  system,  such  that  the  transformations  in 

G,  when  extended  to  complex  values  of  the  variables  (x,t,u),  are  analytic.  Let  G 

o 

denote  the  projected  group  action  on  (x,t) -space.  Assume  further  that  the  action  of  GQ 
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on  some  domain  DQ  is  regular  in  the  sense  of  Palais,  so  that  all  the  G-invariant 
solutions  of  A  =  0  defined  over  are  found  by  integrating  a  system  of  ordinary 

differential  equations#  A/0  =  0,  defined  over  the  image  MQ  of  DQ  in  the  quotient 
manifold  M. 

Theorem  3.4.  Suppose  A  =  0  is  a  fl-completely  integrable  system  of  partial 
differential  equations  in  the  Banach  space  8  with  an  analytic,  regular,  projectable,  one- 
parameter  symmetry  group  G.  _If_  u  =  f(x,t)  is  a  G-invariant  solution  of  A  =  0  with 
initial  data  lying  in  8/  then  the  combination  corresponding  to  Q  of  the  solution  of  the 
reduced  system  of  ordinary  differential  equations  is  meromorphic  in  MQ ,  the  image  of 


C  x  Q  _in_  M. 

Proof.  Since  GQ  is  analytic,  the  orbits  of  0Q  in  the  (x,t) -plane  must  be 
analytic  curves.  If  the  solution  of  the  reduced  equation  had  a  singularity  other  than  a 
pole  on  Mq,  the  corresponding  G-invariant  solution  vrould  have  a  similar  singularity  along 
the  orbit  corresponding  to  the  singular  point.  This,  however,  would  contradict  Theorem 

2.6. 

Thus  Theorem  3.4,  in  a  certain  restricted  sense,  states  that  the  reduced  equation  for 
the  G-invariant  solutions  must  be  of  Painleve  type.  However,  since  the  initial  data  for 
the  G-invariant  solutions  must  lie  in  8 ,  it  is  not  for  every  solution  of  the  reduced 
equation  that  Q  is  required  to  have  only  poles  for  singularities.  In  effect  we  can 
consider  only  those  solutions  which  either  decay  sufficiently  rapidly  at  +  00  along  the 
real  axis,  or  are  periodic  along  the  real  axis.  This  restriction  seems  inescapable  given 
the  particular  method  of  proof.  It  would  be  extremely  interesting  to  remove  these 
restrictions  and  prove  the  conjecture  of  the  introduction  in  fully  generality. 


4.  Applications 

I  *  The  Generalized  KdV  Equations 
Consider  the  equation 

u  +  uPu  +  u  *  0  ,  (4.1) 

t  X  XXX  ' 

where  - p  is  a  nonnegative  integer.  In  Example  3.3,  the  equation  governing  the  scale- 
invariant  solutions  was  derived.  However,  this  third  order  ordinary  differential  equation 
is  rather  complex  to  analyze  in  full,  and  we  therefore  apply  our  results  to  a  simpler  class 
of  self-similar  solutions,  namely  the  travelling  wave  solutions.  Here  the  symmetry  group 
is 


G  :  ( x , t ,u )  ♦  (x  +  cX,t  +  X,u),  X  t  R  , 

c 

where  c  denotes  the  velocity  of  the  wave.  The  invariants  of  Gc  are  £  *  x  -  ct,u,  and 
the  reduced  equation  for  Gc-invariant  solutions  takes  the  form 

u ' • •  +  uPu 1  -  cu '  =0  , 

primes  denoting  derivatives  with  respect  to  This  can  be  integrated  once: 

u* 1  *  — uP+1  +  cu  +  \  d  . 

p  +  1  2 

Multiplying  by  u' ,  a  further  integration  yields 


<u')2  = 


-2 


p+2  2 

u  +  cu  +  du  +  e  , 


(4.2) 


(p  +  D(p  +  2) 

for  some  constants  d,e.  Thus  the  general  travelling  wave  solution  will  be  expressed  in 
terms  of  the  hyperelliptic  function  corresponding  to  the  square  root  of  the  (p  +  2)-th 
order  polynomial  on  the  right  of  (4.2).  The  following  two  results  characterize  the 
singularities  of  the  solutions  of  (4.2). 

Theorem  4. 1  (Painleve's  Theorem).  Consider  the  ordinary  differential  equation 

G(u’ ,u,£)  =  0  , 

where  G  is  a  polynomial  in  u1  and  u,  and  analytic  in  £.  Then  the  movable 
singularities  of  the  solutions  are  poles  and/or  algebraic  branch-points. 

Theorem  4.2.  Consider  the  equation 

(u' )2  -  R(u)  ,  (4.3) 


where  R  is  a  rational  function  of  u.  Then  the  solutions  of  (4.3)  are  all  meromorphic  in 
C  if  and  only  if  R  is  a  polynomial  of  degree  not  exceeding  4. 
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The  proofs  may  be  found  in  Ince,  [23],  and  Hilie,  [72,  p.  6B3] .  Note  that  if  u  has 

an  algebraic  branch  point,  so  also  does  any  linear  combination  of  u  and  its 

derivatives.  Therefore,  for  (4.1)  to  be  linearly  completely  intearable,  (4.2)  must  satisfy 
Theorem  4.2.  Thus  p  -  0,1,  or  2,  and  in  these  cases  the  solutions  are  given  by  elliptic 

or  trigonometric  functions.  Note  that  o  =  0  corresponds  to  the  linear  case,  p  =  1  to 

the  KdV  equation,  and  p  =  2  to  the  modified  KdV  eauation,  all  of  which  are  known  to  hp 
intearable  by  inverse  scattering. 

To  complete  the  demonstration  that  the  generalized  KdV  equations  are  rot  linearly 
completely  intearable  for  p  *  0,1,2,  we  must  place  the  complete  integrability  in  a 
suitable  Banach  space  8,  and  to  do  so  we  check  the  asymptotic  behavior  of  the  travel lina 
wave  solutions  at  +  °°.  I*  we  reauire  that  u'ux  *  0  as  |x[  °°,  then  d  -  e  -  0  in 

(4.2).  Moreover  the  polynomial  on  the  right  of  (4.2)  now  has  a  double  zero  at  u  =  0  and 
a  simple  zero  at  u^  =  [—  ( p  +  1 )  ( p  +  2)c]^^ .  Standard  technicrues,  cf,  [47],  allow  us  to 
conclude  the  existence  of  travelling  wave  solutions  with  positive  velocities  decaying 
exponentially  for  ! x !  ♦  »,  and  reaching  an  extreme  value  of  uQ.  Thus  for  p  odd,  the 
travelling  waves  are  humps  with  uq  the  peak  value,  while  for  p  even,  both  humps  and 
trouahs  occur.  The  important  point,  however,  is  the  exponential  decay  of  these  waves  for 
|x!  •+  w,  and  the  fact  that  for  p  *  0,1,2,  thev  have  complex  nonpolar  singularities.  If 
therefore  we  take  for  the  Banach  space  8  a  snace  of  functions  vanishing  exponentially,  we 
have  shown  that  the  generalized  KdV  equations  are  not  linearly  completely  integrable 
in  8  for  n  *  0,1,2,  and  this  completes  the  demonstration  that  these  eauations  can  be 
solved  bv  inverse  scattering  only  when  p  =  0, 1  or  2.  This  result  is  in  accordance  with 
numerical  evidence,  [15],  that  only  in  these  special  cases  do  the  eauations  have  soliton 
solutions . 

II .  Nonlinear  Kl  e  in -Gordon  Aquations 

Consider  the  nonlinear  V1 e in -Cordon  eauation  in  characteristic  coordinates 

u  =  f* (u)  ,  (4.4) 

xt 

where  f  is  an  analytic  function  of  u,  real  for  real  u,  and  prime  denotes 

■lor  i vat  i  ^ 0  cases  we  wil.1  h<=  nost  interested  in  are  when  f  is  a  polynomial  or  a 


-22- 


finite  sum  of  exponential  functions.  We  will  determine  necessarv  conditions  on  f  for 


(4.4)  to  be  linearly  completely  integrable  by  analysis  of  tbe  singularities  of  the 

travelling  wave  solutions.  If  c  is  tbe  velocity#  £  =  x  -  ct,  then  the  reduced  equation 

for  the  G  -invariant  solutions  of  (4.4)  is 
c 

-cu' 1  =  f* (u)  .  (4.5) 

Multiplying  (4.5)  by  u'  and  intearating  yields 

-  |  (u* )2  =  f(u)  +  k  (4.6) 

for  some  constant  k.  For  simplicity  we  shall  assume  that  k  can  K  chosen  so  that  u^ 

(real)  is  a  simple  or  double  zero  of  f(u)  +  k  and  there  is  a  second  consecutive  simple  or 

double  zero  for  some  real  u2«  This  assumption  ensures  that  the  initial  data  u(x,0)  can 
be  chosen  to  lie  in  a  suitable  Banach  space  8: 

i)  if  u1  and  u2  are  simple  zeros#  so  that  a  solution  of  (4.6)  oscillates 

between  u^  and  u^#  we  take  8  to  be  a  space  of  periodic  functions; 

ii)  if  u.j  is  a  double  and  u2  a  simple  zero#  so  that  a  solution  of  (4.6)  decays 
exponentially  to  u^  as  |£|  ♦  °°#  we  take  8  to  be  a  space  of  functions  exponentially 
converging; 

iii)  if  and  u2  are  double  zeros,  so  that  a  solution  of  (4.6)  tends 

exponentially  to  u^  as  £  ♦  00  and  to  u2  as  £  ♦  (or  vice  versa),  we  can  again 

take  8  to  be  a  space  of  functions  exponentially  converging,  but  to  different  limits. 

The  following  theorem  (stated  in  the  context  of  (4.4)  although  it  applies  generally) 
is  an  immediate  conseauence  of  considering  a  linear  combination  of  u  and  its 
derivatives.  It  tells  us  what  si naularities  are  possible  for  solutions  of  linearly 
completely  integrable  equations. 

Theorem  4.3.  Suppose  for  some  constant  k  that  tbe  analytic  function  f(u)  +  k  has 
two  consecutive  simple  and/or  double  zeros  on  the  real  axis.  Then#  if  the  nonlinear  Klein- 
Gordon  equation  (4.4)  is  linearly  completely  integrable  in  the  relevant  Banach  space 
indicated  above#  it  must  be  the  case  that  any  solution  of  (4.6)  (with  c  having  the 
opposite  sign  to  f(u)  +  k  between  the  zeros)  has  as  singularities  only  poles  or 
logarithmic  branch-points • 


A  logarithmic  branch-point  is  by  definition  a  singularity  such  that  some  linear 
combination  of  derivatives  has  a  pole.  It  arises  in  practice  if  the  scattering  operator 
L  depends  only  on  ux'uXx'***'  so  Qlu]  in  turn  depends  only  on  derivatives,  and  to 

demonstrate  that  this  situation  can  indeed  rise,  consider  the  sine-Gordon  equation 

uxt  =  sin  u  . 

It  was  indicated  in  Example  2.3  that  this  is  completely  integrable,  and  to  examine  it  in 
the  context  of  Theorem  4.3  we  take 

f(u)  =  -cos  u,  k  =  0  . 

The  solution  of  (4.6)  is  then 

/2  sin(^  u)  «  sn{c  ^2  (£  +  6)}  , 

where  sn  is  the  Jacobi  elliptic  function  with  modulus  k  *  1//2  ,  cf.  [9],  This  is  well 
defined  for  c  >  0.  Now  sn  has  simple  poles  on  a  certain  rectangular  lattice  in  C,  and 
so  u  has  logarithmic  singularities  at  these  lattice  points.  The  reason  for  the 
appearance  of  these  nonpolar  singularities  is  the  fact  that  ux  rather  than  u  appears  in 
the  scattering  operator  L,  We  note  that  ux  on  the  other  hand  does  have  only  poles  for 
singularities. 

Theorem  4.4.  Suppose  that  f(u)  is  a  rational  function,  real  for  real  u  and  such 
that,  for  some  k,f(u)  +  k  has  two  consecutive  simple  and/or  double  zeros  on  the  real 
axis.  If  the  Klein-Gordon  equation  uxt  =  f’(u)  is  linearly  completely  integrable,  then 
f  is  a  polynomial  of  degree  not  exceeding  4 . 

The  proof  is  immediate  from  Theorems  4.1-2. 

To  discuss  the  case  where  f  is  a  polynomial  of  degree  <  4,  one  can  try  other 
similarity  solutions  of  (4.4),  or  else  quite  different  tests.  For  example,  it  can  be 
shown,  (12],  that  when  f  is  of  degree  >2,  so  that  f’  is  nonlinear,  (4.4)  has  only 
finitely  many  polynomial  conservation  laws,  while  a  theorem  of  Gel’fand  and  Dikii,  [17), 
[18],  states  that  if  a  system  of  partial  differential  equations  has  a  Lax  representation, 
then  there  are  an  infinite  number  of  polynomial  conservation  laws. 
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Next  we  consider  the  case  where  f  is  a  finite  sum  of  exponential  functions 


m  a  .u 

f{u)=  l  c.e  3  ,  c.  ,a.  e  C  . 
j-0  3  33 

For  simplicity,  we  restrict  our  attention  to  the  case  where  a_.  =  n^a  for  some  a  e  C  and 
some  rational  numbers  n j .  By  dividing  a  by  the  common  denominator  of  the  n^,  we  may 
assume  the  n^  are  integers*  Now  let  v  ■  exp(au),  so  that  v*  -  avu*  .  Thus  v 
satisf ies 


- (v*  )2  =  Z  c .  v  3  .  (4.7) 

2a  3 

Note  that  Theorem  4.2  cannot  be  applied  here  since  v  may  have  singularities  not  shared 
by  u.  However,  since  u'  =v’/o iv,  it  is  necessary  to  find  conditions  on  (4.7)  such  that 
the  function  v'/v,  for  solutions  v,  has  no  movable  algebraic  branch-points.  This 
requires  a  more  detailed  investigation  of  the  proof  of  Theorem  4.2.  It  suffices  for  our 
purposes  to  note  the  following: 

Lemma  4.5.  Consider  the  ordinary  differential  equation 

<v'  )2  =  v-nP(v)  , 

where  P  is  a  polynomial  with  P(0)  #  0  and  n  is  a  positive  integer.  Then  for  any 

c  c  there  is  a  solution  v  with  algebraic  branch-point  at  This  solution  has  a 

Puiseux  expansion 


00 


with  a1  *  0,  and  the  rational  number  r  is  given  by 
i)  r-(m+1)”1  if_  n  =  2m, 
ii)  r  -  2(2m  +  3)  ^  if  n  =  2m  +  1  . 

The  proof  of  this  result  can  be  inferred  from  Hille,  [22,  pp.  681-682). 
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Lemma  4.6 .  Suppose  v  has  an  algebraic  branch-point  at  5  •  Then  v '  / v  has  no 

branch-point  at  t.  if  anti  only  if  v(  C)  -  ( f.  -  t  )rf(K )  for  r  rational  and  f 

o  ■  ■  o  -  ■  — -  ■  ■ 

meromorphic  at  £  « 

Proo*.  Assume  without  loss  of  qenerality  that  £  =0,  Let  v  have  the  Puiseux 

-  o 

expansion 


v (  £) 


r^r 


V 


where  n  is  an  inteqer  and  a  *  0.  Let  a  he  the  first  nonzero  coefficient  for  which 

o  v 

kr  is  not  an  inteqer,  if  such  exists.  Now 


v 


b.c3r 


j-0 


where  b  -  a  and  the  first  nonzero  coefficient  h.  with  ir  not  an  inteqer  is 
o  o  j 


«2 

b  =  -a  a 
k  k  o 


Furthermore 


v' 


jr.r- 1 


jr 


Therefore 


iL 

V 


r3  r 


j-0 


and  the  coefficient  of  £ '  is 

c,  *  -nra  a  1  +  (m  +  k)ra  a"1  , 
k  k  o  Vq 

whjch  vanishes  only  when  a^  *  0.  This  proves  the  lemma. 


Proposition  4.7.  Consider  the  ordinary  differential  equation 


Given  £  £  C,  there  exists  a  solution  v  of  (4.8)  such  that  v»/v  has  an  algebraic 

branch-point  at  £  ,  unless  (4.P)  is  of  the  special  form 


( v* 


C  .  V 

3 


jV  +  2 


(4.9) 


for  some  integer  k . 

Proof.  Let  £  =  n  and  assume  b.  *  0,  b  *  0.  By  Lemma  4.«>  all  solutions  must  be 

'  o  N  — n 

j* 

of  the  form  v(  £)  *  £  f(£)  with  r  rational  and  f  meromorphic  at  0  if  we  are  to  avoid 

an  algebraic  branch-point  for  v’/v.  Thus 

2  2r  -1  2 

(V)  =  £  ( r  £  f  +  fM  , 

and 

vj  =  r,jrfj  , 

so  that,  eauating  the  fractional  powers  of  £,  we  see  that  b_.  *  0  unless  jr  *  2r  +  i 
for  some  intener  i.  If  n  >  0,  it  follows  from  Lemma  4.5  that  b_.  =  0  unless 
i )  j  =  2  nod  ( m  +  1 )  for  n  =  2m  , 
or  ii)  2j  E  4  nod  (2m  +  3)  for  n  =  2m  +  1  . 

In  particular,  the  only  neaative  values  of  j  which  satisfy  these  congruences  are 
1  -V^n  and  -n,  the  first  value  occurring  only  when  n  is  even. 

Next  set  w  *  1/v.  Then  (4.8)  becomes 


( w* ) ^  «  :  b,w^  ^  • 

j-n  1 

Since  w*/w  =  -v'/v,  w  must  satisfy  the  same  conditions  as  v.  Therefore,  if 

M  >  4,  b.  =  n  unless 
3 

i )  j  i  2  mod  tv  -  1 )  if  v  =  21' , 

or  ii  )  2j  "  4  no'4  <  2*'  -  2 )  if  *’  =  ?'*  +  1  . 

The  only  positive  values  of  i  satisfyino  these  are  N,  V2N  +  1  and  2,  the  second  only 
if  *:  is  ef'en .  'Y>nnar  i son  of  the  two  sets  of  congruences  then  shows  that  (4.8)  must  be  of 
the  required  'om. 
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Theorem  4»8.  Suppose  f(u)  is  a  linear  combination  of  exponential  functions  e  J 
with  *  rijCt,  rational/  a  complex.  Suppose  further  that  f(u)  is  real  for  u 
real,  and  that,  for  some  real  X,  f(u)  +  k  has  two  consecutive  simple  and/or  double  zeros 
on  the  real  axis.  If  the  Klein-Gordon  equation  uxt  -  f'(u)  is  linearly  completely 
integrable/  then  f  must  be  of  the  special  form 

f(u)  =  [  c.ei&U  ,  (4.10) 

3--2  3 

where  0  is  a  rational  multiple  of  a. 

It  is  interesting  that  the  form  (4.10)  for  f  includes  the  double  sine-Gordon 
equation 

u  *  a  sinotu  +  b  sin(-J-  au)  , 
xt  2 

for  which  numerical  studies  of  Dodd  and  Bullough,  (11]#  indicate  the  existence  of  soliton 
solutions.  A  recent  preprint  by  Gibbons  and  Fordy,  [20],  contains  the  result  that  the 
special  case  of  (4.10)  when  f(u)  -  e2u  +  e“u  does  have  a  Lax  representation,  but  it  is 
not  known  whether  the  result  extends  to  a  general  function  f(u)  of  the  form  (4.10). 

III.  Model  Wave  Equations  of  Whitham  and  Benjamin 
The  integro-dif ferential  equation 

ut  +  uux  +  tf[ux)  *  0  ,  (4.11) 

where  H  is  the  integral  operator 

00 

H[f ]  (x)  -  /  HU  -  y )  f  ( y )  dy  , 


was  proposed  by  Whitham,  [46J,  [47],  as  an  alternative  to  the  KdV  equation  for  long  waves 
in  shallow  water  which  could  also  model  breaking  and  peaking.  Here  H  is  taken  to  be  the 
Fourier  transform  of  the  desired  phase  velocity  c{k),  where  k  is  the  wave  number.  Of 
particular  interest  is  the  case 

c(k>  -  -r-2 — r  ,  v  >  0  , 

v  +  k 
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•  v 


so  that 


ho.)  -  f  e-v,xl  . 

2  v 

2  2 

Note  that  H  is  the  Green's  function  of  the  operator  D  -  v  *  V  so  that  (4.11)  is 
equivalent  to  the  differential  equation 

P(u+uu]+u*0.  (4.12) 

t  X  X 

It  can  be  shown,  [15],  that  (4.12)  possesses  travelling  wave  solutions  u,  with  |u|  +  o 
as  |x|  *►  «,  and  amplitudes  between  0  and  some  maximum  height.  Computer  studies 
indicate  that  these  waves  may  be  solitons,  i.e.  interact  cleanly.  One  possibly  undesirable 
feature  of  (4.11)  is  the  extremely  fast  propagation  of  short-wave  components,  and  for  this 
reason  Benjamin,  Bona  and  Mahony,  [7],  proposed  the  alternative  model 

(4.13) 

Again,  in  the  special  case,  (4.13)  can  be  rewritten  as 

(4.14) 

In  general,  we  will  let  V  be  any  constant  coefficient  linear  differential  operator 


u  +  uu  -  H[u  1  -  0  . 

t  X  t 


p[ut  +  uu^]  -  u  *  0  . 


n 

=  I 

i*n 


c.  D1 ,  c  /  0  . 
i  n 


We  show  here  that  the  model  equations  (4.12),  (4.14)  cannot  be  integrable  by  inverse 
scattering  methods.  As  usual,  consider  the  travelling  wave  solutions  of  these  ecruations. 

If  c  denotes  the  velocity,  then  the  reduced  ecruation,  after  integration,  is 

Plj  (u  -  c) 2 )  +  a(u  +  d)  -  0  .  <4.1*1 

Here  d  is  a  constant  of  integration,  a  *  1  in  the  Whitham  model,  a  =*  c  in  the 
Benjamin  model,  and  D  now  denotes  d/d£,  £  *  x  -  ct.  since  n-th  order  eauations 
of  Painleve  type  have  not  been  classified,  we  resort  to  Painleve’s  original  "OMnethod"  to 
analyze  the  singularities  of  the  solutions  of  (4.15).  The  basic  result  is  found  in  Ince, 
(23,  p.  319). 

Lemma  4.9.  Suppose  MutC,,a)  *  0  is  an  analytically  parametrized  family  of  ordinary 
differential  equations  for  a  in  some  domain  C:  conta ining  0  as  an  interior  point.  If 
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the  general  solution  u<  r. ,  a)  is  uniform  in  £  for  ct  c  P  {0},  then  it  will  be  uniform 
for  a  =  0. 

In  our  case,  let  r,  ~  +  a£.  Then  if  we  consider  u  as  a  function  of  (4.  IS) 

becomes 

(c  Pn  +  or  1  +  ...  +  oV  )  r-i  (u  -  c)2]  +  ana(u  +  d)  »  0  , 

n  n- 1  o  7. 

where  D  now  denotes  d/d£.  For  ct  =  0,  this  reduces  to 

n  1  2 

D  (—  (u  -  c)  )  =  0  , 

the  solution  of  which  is 

u  =  c  +  /P  ( O 
n 

for  an  arbitrary  polynomial  Pn  of  degree  <  n  -  1.  This,  for  appropriate  Pn ,  has  an 

algebraic  branch-point  at  C  =  0,  so  that,  by  the  lemma,  solutions  of  (4.15)  must  also 

have  non-loaarithmic  branch-points.  'This  involves  a  slight  extension  of  the  lemma  above, 

but  it  is  easy  to  infer  its  truth  from  the  proof  given  by  Ince . )  If  these  solutions  also 

satisfy  decay  or  periodicity  properties,  Theorem  3.4  (together  with  Theorem  4.3)  shows  that 

model  equations  (4.12),  (4,14)  cannot  be  linearly  completely  integrable.  In  particular, 

2  2 

Whitham's  equation  with  V  -  D  -  v  is  not  integrable  by  inverse  scattering. 

I v •  The  PFM  Equation 
The  ecruation 


u.  +  uu  -  u  ^  =  0  (4.16) 

t  X  xxt 

was  proposed  by  Benjamin,  Pona  and  Mahonv,  [7? ,  as  an  alternative  model  to  the  KdV  eauation 
for  the  description  of  long  waves  in  shallow  water.  In  [36]  it  was  shown  to  possess  only 
three  independent  conservation  laws,  and  therefore  by  the  results  of  Gel'fand  and  Dikii 
cannot  be  completely  integrable  Our  consideration  of  this  example  runs  into  difficulties 
because  the  self-similar  solutions  do  not  satisfy  any  decay  or  periodicity  properties,  and 
the  functions  9  we  can  allow  are  limited,  but  we  will  indicate  the  method  here. 

First  we  note  that  (4.16)  admits  the  symmetry  group 


f  > 


x , 


v  ,  -XX, 

O )  m  ( x  ,e  t  ,e  u)  ,  X  c  R  , 


jf  scale  transformation t,  !- variants  o*  0  are  provided  by  x  and  w  =  tu,  for  t  >  0, 

*or  "i-i'.-MMant  solutions  is  then 


sr1  v 


ice*:  evja 
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w' *  +  ww'  -  w  =  0  ,  (4.17) 

the  primes  denoting  derivatives  with  respect  to  x.  It  can  be  readily  checked,  by  the 
procedure  in  Ince,  [23],  that  (4.17)  is  not  of  Painleve  type.  Indeed,  it  is  of  type  i(b) 
on  page  330  of  Ince.  Applying  the  a-method  as  Ince  does,  one  can  readily  check  that 
branch-points  appear,  although  possibly  only  logarithmic,  and  this,  granted  the  existence 
of  a  suitable  Banach  space  8,  would  show  that  the  BBM  equation  is  not  Q-completely 
integrable  for  Q,  say,  the  identity. 

However,  a  closer  investigation  of  the  behavior  of  the  real  solutions  of  (4.17)  is 
required.  Since  x  does  not  appear,  it  can  be  integrated  to  yield 


w*  -  Vo  w2 

( 1  -  w* )e  *  ce  2  .  (4.18) 

In  principle,  this  equation  can  again  be  integrated  by  solving  for  w'  in  terms  of  w.  To 
investigate  the  solutions  qualitatively,  note  that  w'  =  0  if  and  only  if  w  =  2  log  c, 
c  >  1.  The  only  double  root  is  when  c  **  1 ,  and  only  in  this  case  do  solutions  decay  at 
+®  or  «oo.  However,  it  is  readily  seen  that  a  solution  decaying  at  one  endpoint  cannot 
decay  at  the  other,  nor  are  periodic  solutions  possible.  Thus  we  are  unable  to  apply  our 
results  to  this  case. 

V.  Lax  Pairs  of  Composite  Order 

Gel*fand  and  Dikii,  [17],  [18],  succeeded  in  classifying  all  Lax  pairs  of  differential 

operators  of  the  following  special  type.  Let 

Ln  “  D"  +  Un-Z0""2  +  ”•  +  U1D  +  “0 

40 

be  a  scalar  differential  operator  of  order  n  with  u  =  (Uq,*..,u  independent  C 

functions,  and  D  =  d/dx.  They  showed  that  for  each  integer  m  not  a  multiple  of  n, 

there  is  a  differential  operator 

P  =  0m  +  p^  ^  -Dm”2  +  ...  +  I*.  •  D  +  P 
m  Mn,m-2  *m,1  Km,0 

of  order  m,  the  pm  i  biing  polynomials  in  the  u^  and  their  derivatives,  such  that  the 


Lax  representation 


9L 

ar  ■  'VV 
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is  a  genuine,  nontrivial  system  of  evolution  equations 


ut  "  Km(u>  •  (4.19) 

Moreover,  the  Pm  are  unique  if  we  require  the  coefficients  pm  ^  to  have  no  constant 
term. 


Consider  the  stationary  solutions  of  the  system  (4.19),  i.e.  those  in  which  u  is 
independent  of  t.  These  satisfy  the  system  Km(u)  =0,  or  equivalently,  the  "stationary 
Lax  representation" 


tPm.Ln]  =  0  .  (4.20) 

Theorem  4.10.  If  the  orders  n,  m  of  the  operators  Ln,  Pm  in  the  Lax 
representation  of  (4.19)  are  not  relatively  prime  integers,  then  stationary  solutions  of 
(4.19)  with  arbitrary  singularities  in  the  complex  plane  exist. 

Proof .  Let  k  >  1  be  the  greatest  common  divisor  of  m  and  n.  Consider  the 
operator 


Mk  =  Dk  +  vk-2Dk~2  +  +  viD  +  vo  ' 

whose  coefficients  v^(x)  are  sufficiently  differentiable  for  x  e  R  but  are  otherwise 
arbitrary  functions.  Then 


un  ,0 


=  <MV> 


n/k 


^m,0 


(M.  ) 


m/k 


obviously  satisfy  the  stationary  Lax  representation  (4.20)  and,  moreover,  using  the 
formalism  of  Gel* f and  and  Dikii,  it  is  easy  to  prove  that  pm  ^  is  derivable  from  Ln  ^ 
via  the  same  formulae  as  gave  Pm  from  Ln.  Therefore  each  such  gives  a  stationary 

solution  of  the  evolutionary  system  (4.19). 

Now  suppose  that  Ln  is  any  such  operator,  where  n  is  a  composite  number.  If  there 
exists  a  Gel  * fand-Levitan  type  of  integral  equation  for  solving  the  inverse  problem  for  the 
operator  Ln,  then  Theorem  3.4  would  imply  the  meromorphic  character  of  the  group- 
invariant  solutions  of  the  evolutionary  system  (4.19),  using  similar  arguments  to  those 
used  in  the  integration  of  the  Korteweg-de  Vries  equation.  This,  however,  is  in 
contradiction  to  Theorem  4.10  for  the  case  of  time-invariant  solutions.  (The  relevant 
symmetry  group  is  just  translation  in  t.)  This  indicates  that  such  a  differential 
operator  of  composite  order  does  not  have  an  inverse-scattering  formalism  in  the  sense  that 
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the  Schrodinger  operator  does  -  either  no  such  Gel  * fand-Levitan  equation  exists,  or  the 
assumptions  regarding  analyticity  are  not  justified.  Indeed,  we  know  of  no  such  Gel' fand- 
Levitan  equation  for  any  operator  of  composite  order,  e.g.  for  order  n  *  4. 
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A  completely  integrable  partial  differential  equation  is  one  which  has  a 
Lax  representation,  or,  more  precisely,  can  be  solved  via  a  linear  integral 
equation  of  Gel ' f and-Levitan  type,  the  classic  example  being  the  Korteweg- 
de  Vries  equation.  An  ordinary  differential  equation  is  of  Painleve  type  if 
the  only  singularities  of  its  solutions  in  the  complex  plane  are  poles.  It 
is  shown  that,  under  certain  restrictions,  if  G  is  an  analytic,  regular 
symmetry  group  of  a  completely  integrable  partial  differential  equation,  then 
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20.  ABSTRACT  -  Cont'd. 


the  reduced  ordinary  differential  equation  for  the  G-invariant  solutions  is 
necessarily  of  Painleve  type.  This  gives  a  useful  necessary  condition  for 
complete  integrability ,  which  is  applied  to  investigate  the  integrabili ty  of 
certain  generalizations  of  the  Korteweg-de  Vries  equation,  Klein-Gordon  equa¬ 
tions,  some  model  nonlinear  wave  equations  of  Whitham  and  Benjamin,  and  the 
BBM  equation. 


